On the unlocking of color and flavor in color-superconducting quark matter 
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The role of the strange quark mass for the phase structure of QCD at non- vanishing densities is 
studied by employing a recently developed self-consistent truncation scheme for the Dyson-Schwinger 
equations of the quark propagators in Landau gauge. Hereby the medium modification of the 
effective quark interaction by the polarization of gluons is implemented. Taking into account this 
effect results in significantly smaller dynamical quark masses at the Fermi surface. Due to this 

\£l ' reduction the color-flavor locked phase is always the preferred color-superconducting phase at zero 

f"^ ' temperature and for a realistic strange quark mass. 

o ' 
<N 

Or 

! The possible occurrence of different phases in Quantum Chromo Dynamics (QCD) at sufficiently high densities 
and small temperatures has attracted a lot of interest in recent years (for corresponding reviews see e.g. refs. 0,01)- 
Within models, and depending on various parameters and constraints, a surprisingly rich structure of different color- 
superconducting phases has been studied. The pairing pattern of quarks strongly depends on the relative position of 
the Fermi surfaces of the different quark species. Therefore the dynamically generated constituent masses as well as 
the effective chemical potentials are the relevant quantities for a study of color-flavor unlocking. 

The color-flavor locked (CFL) phase is expected to be the ground state for a three flavor color-superconductor at 
asymptotically large densities |3| ■ This raises the question of what happens to this state when the density is lowered to 
' more realistic values. Since the CFL phase is neutral for degenerate quarks this depends mainly on the mass splitting 
between the strange quark and the two light (up and down) quarks. The unlocking of color and flavor has been 
^sO • investigated within models of the Nambu & Jona-Lasinio (NJL) type by varying the strange quark mass It has 

been concluded that at some densities the CFL phase is disfavored compared to the 2-flavor color superconductor (2SC) 
phase in a self-consistent treatment 0. The region of the 2SC phase possesses a rich (sub-)structure if furthermore 
fS |' neutrality conditions and /3-equilibrium are taken into account [7fl. The CFL phase and the boundary for unlocking 
JL , are, however, only mildly affected. 

These results, obtained in NJL-type models, are rather sensitive to model parameters. For more reliable statements 
a non-perturbative approach, which is directly based on the QCD degrees of freedom and recovers the correct vacuum 
1 results as well as the proper weak coupling behavior at asymptotically large densities, is needed. 

For the description of QCD vacuum properties, truncated Dyson-Schwinger equations (DSEs) in Landau gauge 
have been successfully applied (for corresponding reviews see refs. 0). They naturally describe gluon [j]} as well as 
quark |10| confinement. They have furthermore been extended to non- vanishing chemical potentials with the correct 
behavior in the weak coupling limit Therefore this method will be employed here for the investigation of the 

color-flavor unlocking. Within this framework, we present self-consistent solutions of the quark propagators in the 2SC 
and CFL phase for realistic values of the strange-quark current mass. We include energy- and momentum-dependent 
dressing functions for the most general parameterization in both Dirac and color-flavor space. For the reasons stated 
above, we neglect for the time being the neutrality conditions in the present investigation of color-flavor unlocking. 
In contrast to studies employing NJL-type models, we have also implemented a medium modification of the effective 
quark interaction taking into account the medium contribution to the polarization tensor which enters the gluon 
propagator. It has to be noted that a proper description of the weak coupling limit requires the inclusion of this 
effect. The vacuum part of the effective quark interaction is taken either from a solution of coupl ed g luon, ghost and 
quark propagator DSEs or from a fit to lattice data for the gluon and the quark propagators jl2lll3| . 

In the chiral limit and at moderate values of the quark density, the gap functions at the Fermi surface in such 
a calculation acquire sizeable values 11], comparable to those obtained within NJL-type models. These results are 
found to be rather insensitive to the effective quark interaction, at least within the range of uncertainty of the effective 
coupling. The dynamically generated mass functions, on the other hand, turn out to be significantly smaller than 
in the vacuum. This is due to the medium modification of the interaction (which is usually neglected in NJL-type 
investigations). This effect of smaller dynamical quark masses implies that the Fermi surfaces of the quarks are closer 
to each other than expected and eventually leads to the most important result presented here: the CFL phase is, at 
zero temperature and realistic values for the strange quark current mass, always energetically preferred compared to 
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the 2SC phase. 

This paper is organized as follows: In section [n] we briefly present the employed theoretical framework, leading 
to a fully self-consistent treatment of the CFL phase for a range of strange quark masses. A physical interpretation 
of the parameterization used is provided. In section IIIII we discuss the numerical results including those for Fermi 
momenta, mass and gap functions, etc.. Most importantly, we provide a prediction for the critical value for the strange 
quark current mass as a function of the chemical potential. Finally, we summarize and conclude in section IIVI The 
color-flavor basis employed is given in the Appendix. 



explicitly leads to 



II. THE LANDAU-GAUGE QUARK PROPAGATOR AT NON- VANISHING CHEMICAL POTENTIAL 

A. The truncated quark DSE at non-vanishing chemical potential 

In ref. a closed truncated DSE for the Landau-gauge quark propagator at non- vanishing chemical potential 
has been derived and solved in the chiral limit. In this work we are following the scheme and notations presented 
there. In the Nambu-Gor'kov basis (see ref. for a corresponding review) the normal quark DSE is coupled to an 
equation for the gap functions. To describe their structure we define 

Soip) = ( S °o P) So(p) = -cs+{-p) t c ) ' (1) 

* [p) ~ { T+(p) S-( P ) = -CS+(- P ) T C ) ' {l) 

(^ +{p) ^ ip) ] (3) 
MP) - ( v $+( :p ) S -(p) ) ' {6) 

with So being the bare and S the full Nambu-Gor'kov propagator. £ is the Nambu-Gor'kov self-energy. Sq is diagonal 
in color-flavor space and equal to S^^p) -1 = —i(p4 + in)j4 — ip-j + moj for a quark flavor /. As the Euclidean action 
is real we furthermore have to impose 5> — (p) = 74 < I >+ (p)^74. Introducing appropriate renormalization constants, the 
DSE for the quark propagator 

S-\p) = Z 2 S - 1 (p) + Z 1F n P ) (4) 

r ± = -z 1F (z 2 s^ 1 + z 1F Y^y\ ± s ± , (5) 

S^ 1 = Z 2 S t ~ 1 + Z 1F E ± - Z\ F ^ (ZiS^~ X + 1 $ ± - (6) 

Employing the truncation of ref. [llj], the equations for the self-energy and the gap function then read 

*® = ^Jwr^ s+ ^ (ray + ray) • (7) 

ffi+M z " f diq xtt+m > ( a °^ p £ i a ^)p^ \ ™ 

* w = -zT/1 wr 7 " a {q)lv a [w+m ) ' { ) 

where k = p — q. Here projectors transverse and longitudinal in the medium have been introduced. The functions 
G and F describe the corresponding medium modifications of the gluon propagator, see Eq. (21) of ref. [TT| and are 
calculated, once the coupling is given. Therefore the only input for the quark DSE are the running coupling a s (k ) 
(for its choice see below) and the renormalized quark current masses. 

Before proceeding we note that this quark DSE can also be derived from the thermodynamically consistent effective 
action corresponding to the truncation scheme used |Tlj . Evaluating this action at the stationary point one obtains: 

r[<S] = --Tr P)D>C; / iW - G Ln5 _1 + 2 Tr p,£',c,/,JVG (l - Z 2 S ~ 1 S) + const. 

= ~Tc PiD)Cl fLa. (Z 2 S^ ~ l + Z 1F 5T)) + 

+ jTr p , AcJ (2 - Z2S+S+- 1 - Z 2 S- So - 1 ) + const. (9) 
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Later on we will use this expression to determine e.g. the pressure difference between different phases at a given 
chemical potential. 

Similar to our previous investigation in the chiral limit we employ two different couplings. The running 
coupling determined in DSE studies of the Yang- Mills sector 12], which is labeled in the following as a/(fc 2 ), will 
serve as a lower bound. As is detailed in ref. [12j it underestimates chiral symmetry breaking significantly in the 
abelian approximation, which on the other hand is an appropriate choice if the ground state does not break chiral 
symmetry dynamically or explicitly, see ref. |lt| and references therein. The running coupling extracted from lattice 
QCD data for the quark and gluon propagators 13] will serve as an upper bound, since the medium will probably 
weaken the interaction. It is labeled an (k 2 ). In order to vary the renormalization point in the latter version of the 
coupling, we use the multiplicative renormalizibility of the quark DSE. Note also that the vacuum quark propagator, 
determined within such a scheme |l4j] is in excellent agreement with the corresponding lattice data when the finite 
volume of lattice calculations is accounted for. 

As stated above, taking into account the medium polarization, Debye screening and Landau damping are included. 
Both, screening and damping of the interaction increase for increasing interaction strength. Therefore the generated 
gap and mass (normal self-energy) functions turn out to be much less sensitive against variation of a s (k 2 ) than the 
dynamical mass function in the chirally broken vacuum. Due to this fortunate instance the presented approach has 
considerable predictive power despite the uncertainty in the effective low-energy quark interaction at non-vanishing 
densities. 

The renormalization constants are determined in the (chirally broken) vacuum. Due to the employed vertex con- 
struction, the quark-gluon vertex renormalization constant, Zip, cancels in the resulting renormalized equations. For 
each flavor, we determine the quark wave-function renormalization constant, Z 2 , and the renormalization constant 
Z m , relating the unrenormalized quark mass mo l<2 (A 2 ) at an ultraviolet cutoff A to the renormalized mass m q (v) via 

m , 9 (A 2 ) = Z m (v 2 ,k 2 )m q (v), (10) 

by requiring 

s q(p)\ p 2 =u 2 = -if + m q {v) (11) 

at a renormalization scale v. This corresponds to a momentum-subtraction (MOM) scheme, which results in some- 
what smaller numerical values for the quark current masses at a given renormalization scale (usually taken to be 
2 GeV) . We simply ignore here the difference between MS and MOM masses because the effect is of the order of ten 
percent (when calculated within perturbation theory) and thus much smaller than the uncertainty in the value the 
current masses poj. 

It turns out that, as expected, the mass dependence of the quark wave function renormalization constant, Z 2 , is 
negligible as long as the renormalization scale is much larger than the mass. Therefore, we simply drop this dependence 
and Z 2 is determined once and for all in the chiral limit. To keep the number of parameters as small as possible, we 
work in the chiral isospin limit (i.e. we set the up and down current quark masses to zero) and vary only the strange 
quark current mass. 

In the following we will restrict to isotropic phases. In order to solve the DSE of the quark propagator it is 
advantageous to consider their color-flavor structure first. To get a self-consistent solution we choose suitable sets of 
matrices {Pi} and {Mi} in color-flavor space, such that 

^t(p)Pi, (12) 
(t>tiP)Mi, (13) 

where we have introduced the renormalization-point independent component functions T,f(p) and (f>f(p), which are 
matrix-valued in Dirac space. Full self-consistency is guaranteed in case a basis of all allowed matrices is considered. 
The dimensionality of this basis in a given phase depends on the residual symmetry in color-flavor space. For the 
CFL phase this will be detailed below. 

The Dirac structure of the self-energies in an even-parity phase can be parameterized by |16| 

E+(p) = -^E+ i (p)-^ p E+ i (p) + S+»-n^E+ i (p), (14) 

#(P) = (74M>)+74 4>) + ^>)+M>))75, (15) 

where p = p/\p\, ft = p ■ 7, i/> — w p 74 and ui p — ip^ + /i. We already note that i (jp) turns out to be negligibly small. 
This is to be expected since it has to vanish due to time reversal symmetry at vanishing temperatures in color-flavor 



E+(p) = 


z 2 

Zip 


$+(p) = 


z 2 

ZiF 
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symmetric channels. The gap functions 4>b and 4>d, which vanish in the chiral limit will be discussed in subsection 
111 CI where we provide an interpretation in a simplified setting. 

The system of equations for the self-energies is symmetric under the transformations 

Za/b/cJP4>\p\) s 1/b/c^(-P4,|p1)*, (16) 

££>4>1) - -E+ ,(-p 4 , \P\T, (17) 

which we find to be unbroken in the self-consistent scheme presented here. Therefore, the determinant of the quark 
propagator is positive for p^ — 0, and the Luttinger theorem is applicable for determining the density . 

For the anomalous propagator an inhomogeneous part is missing in Eq.JSJ. Therefore we can choose a global phase 
for the gap functions. This is done such that the equations are invariant under 

<i> + a/d ,mm) ^ <Pt /D ,ii-Pi,\p\r (is) 



04,1$) - -<t>t/c,i(-p*>\p\y, (is) 



and therefore 4>a and <j)£> are real, and 4>b and <pc purely imaginary, for p± = 0. 

Neglecting i: the Dirac structure of the inverse normal propagator can be written as 



S+ \p) = Z 2 S+ t 1 + Z 1F Et = -i?A i (p)-ty p Ci(p)+B i (p), (20) 
which allows to define the mass function for the pairing quasiparticles for a given color-flavor channel q as 

TOn a + St „(p) 

MJp) = ' q , B ' qK "' . (21) 
q ! + 

In the following we will mostly present results for mass functions at p± — 0. To ease notations we will use M q {\p\) — 
MqiPi = 0, \p\) except when stated explicitly otherwise. 



B. The parameterization of the CFL phase 

As its name already expresses, the CFL phase is defined via its symmetry pattern in color-flavor space. For three 
degenerate quark flavors, one considers the SU(3) c+ v symmetry generated by r — A^, with a = 1, ... ,8 and r a , A a 
being the Gell-Mann matrices in flavor and color space, respectively • The quarks therefore are ina3cg)3 = l©8 
representation of this symmetry. In the case of only two degenerate quarks and a strange quark, the symmetry is 
broken down to SU(2) c+ y ® [7(l) c -|_y, which is generated by r a — with a = 1, 2, 3 and 8, i.e. the quarks form a 
1©1©2©2©3 representation. 

For the Nambu-Gor'kov propagator to be invariant under this transformation, we need to require 

U^S + (p)U = S + (p), U T T+(p)U = T + (p), (22) 

for U 6 SU(2) c+ v ® £/(l) c +y. The matrices {Pi} and {Mi}, needed for a self-consistent description are then those 
given explicitely in the Appendix. 

In section lllll we will also present results for the self-energies evaluated at different energies and momenta, in 
particular at the Fermi surface. In contrast to the chiral limit, where the matrices {Pi} can be chosen as constant 
projectors onto irreducible representations of the residual symmetry, the situation is more complex here. In principle, 
every dressing function of the normal propagator and self-energy can be decomposed into irreducible projectors in 
color-flavor space: 

F = F 1 P 1 +F 1 ,Pr+F 2 P 2 + F2P2 + F 3 P 3 , (23) 

where, in particular, P 2 = P7, P 2 = Ps and P3 = Pi — \Pi- However, for every dressing function, the singlets arc not 
protected against mixing and therefore their projectors P\ and Py are in general energy and momentum dependent. 
We therefore find six different Fermi momenta, which are defined as sign changes of the determinant of the propagator, 
det(5(p4 = 0, |p|)), at P4 = and connected to the density by the Luttinger theorem. Three of them correspond 
to the upper 3x3 block-matrix in the ansatz employed (see Appendix) and describe the mixing of the two singlet 
channels as well as the triplet channel. After ordering, these are denoted by p_f,i i: Pf.i 2 and pf.i 3 - In the lower 6x6 
diagonal block- matrix, we find pf,2, Pf,2 and pf,3 corresponding to P7, Ps and Pq, respectively. 
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For the dressing functions of the anomal propagators and the gap functions the situation is similar. With M — 
\ 2 ® t 2 + A5 <g> r 5 + A7 <g> T7 being invertible and fulfilling U T MU — M, we define Mj = MPj. Every dressing function 
can then be written as 

G = GxMi + G V M V + G 2 M 2 + G2M2 + G3M3, (24) 

where M2 = M7, M^ = Ms and A/3 = Mi — 5M2. Furthermore Mi and M x / are again in general energy and 
momentum dependent due to the possible mixing between the singlets. 



C. The gap functions 4>b and 4>d 

Having introduced the most general Dirac structure for the gap functions in an even-parity phase (see Eq.JTSJ) 
their interpretation and in particular their relation to the energy gap in the excitation spectrum is of interest. To this 
end one has to determine the dispersion relations as given by the poles of the propagators. This is a solvable task 
but leads already for degenerate quarks and constant dressing functions to very involved and intricate expressions. 
A physically motivated approximation with a simple interpretation of the gap function will be introduced in this 
subsection and will serve as an illustration. For simplicity, we consider the pairing of only two different quasiparticles 
a and b with inverse propagators in Dirac space, given by = Z 2 S^ a i^ + Z 2 T^, b for vanishing gap functions. 

This is in particular the case for the three lower 2x2 block-matrices in the ansatz presented in the Appendix. 

The gapped propagator S+ in one channel, see Eq. JBJ, is then typically given by 

S+- 1 = S+- 1 - ( f>-S^cP+, (25) 

with cf> + of the form in Eq. (JTSJ. The propagators, being of the form in Eq. (|20|l . can be expressed via the energy- 
projectors for i = a,b,c 



^(-wCi-eEt)^, (26) 



e=± 



where 



Af = \ (i± Wi4 + <*m)), (27) 



Pa ~ <t>% ~ <t>c + 4>l 



and Ei = y/jPAf + B'f, on = Bi/Ei and fa = \p\Ai/Ei. Neglec ting the ^-dependence in the dressing-functions and 
absorbing C a and by rescaling E a , Eb and fa correspondingly J21|, we assume that no pairing between quasiparticles 
and anti-quasiparticles takes place. This reduces Eq. H25fl to 

(-iwC c -B c )7 4 A+ w (-no - E a ) 74 A+ - 0- , 1 A fc '740 + , (28) 

— iu* + Eb 

which is exact in the chiral limit. Note, that it can also be justified, at least in leading order, in a fa 2 //J? expansion. 
Evaluating the roots of the right-hand side, we find the dispersion relation 

E a — Eb ( ( E a + Eb 
- IP4 = ___ ± ^J____ M 

+2ifaj {fafaA - a 2 i(t> B ) - 20d {u 2 faA + foifas) J , (29) 

where we have made use of the fact that 4>a and fan are real and fag and faj are imaginary at p<t = (see the discussion 
in the last subsection). Noting that af + Pf = 1, the above expression suggests to introduce gap functions, which are 
transformed accordingly: 

4>A = Ma - a 2 ifa3, i<f>B — oi 2 fa\ + foifas- (30) 

This transformation depends on the mass parameter of the quasiparticle, which makes evident that the existence of 
the functions fas and fao allows different energy gaps for quasiparticles with different masses. In terms of the 'rotated' 
gap functions fa\ and 4>b we obtain 

~*P4 = ^^±((^±^-^,r + (4>A+^faJ) 2 + (fao~^4>B) 2 y 2 ■ (31) 

In analogy to the situation in the chiral limit, we arrive now at a simple interpretation: One has a chirally symmetric 
pairing via fa\ + i<f>C and a chirality breaking pairing through fao — i<pB- For anti-quasiparticle pairing fa\ and fas 
simply change sign. In the chiral limit we then find, as expected, (j)^ = fa\ ± i4>c and fas = fao = 0. 
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III. RESULTS 

As described above, the quark propagator is highly non-trivial. The Dirac structure of the self-energies and gap 
functions is given by four self-consistently determined functions, respectively, which are in addition functions of \p\ 
and P4. In accordance with the aim of this paper, namely to investigate the phenomenological importance of the 
strange quark current mass, we mainly restrict ourself to the presentation of results at P4 = and some values of 
\p\. This will also turn out to be sufficient to demonstrate the important differences of our results as compared to 
corresponding ones obtained in NJL-type models. As already mentioned, we include self-energy effects, which have 
been analyzed so far only in the weak coupling regime |17|. Note that in the chiral limit the dependence on \p\ has 
already been discussed in ref. |Tl) and the role of a non-trivial ^-dependence in ref. |18| . 

The following results will be presented for the couplings discussed above, ai(k 2 ) and au(k 2 ). As explained, we 
consider them as the limiting cases, which are allowed by the uncertainty within investigations of infrared QCD. We 
will see that gap functions and Fermi momenta are quite insensitive to the coupling used. 

A. Fermi momenta 




Figure 1: Fermi momenta in different channels (see text) at fi = 400 MeV for the coupling ai(k 2 ) (left) and au(k 2 ) (right). 

In Fig.^the results for the Fermi momenta at a chemical potential of fi = 400 MeV as a function of the renormalized 
strange quark current mass m s {v) at a renormalization scale v = 2GeV are presented. For the CFL phase these are 
of course only plotted below the critical value of the strange quark current mass. (For the definitions of the different 
components see section lIlIj| ) Above this critical value the 2SC phase is energetically preferred, and the three different 
Fermi momenta of the 2SC phase are shown: For the gapped red and green, up- and down-quarks and the ungapped 
blue up- and down-quarks, which are both independent of the strange quark mass and for the decoupled strange 
quarks in the unbroken phase (which is also displayed below the critical of m s ). 

In the CFL phase the Fermi momenta are not monotonous functions of m s . This is due to the behavior of the 
vector self-energies and not visible in the mass functions, cf. Fig. [21 

B. Mass functions 

As a function of m s , we display in Fig. [21 the results for the constituent quark mass function in the vacuum at 
vanishing momenta M s ^roken vacuum (0), as well as the quark mass function in the unbroken phase M SiUn b ro k en (PF) 
and the doublet channel in the CFL phase M s ^cfl-i{pf)i both at their respective Fermi momenta and for a chemical 
potential of fi = 400 MeV. Furthermore, the quark mass function in the unbroken phase M s ^ un b ro i- en (2 GeV) for the 
same chemical potential at the renormalization scale v = 2 GeV is given. 

Although the constituent quark mass in the vacuum M s ^ ro k envacuum {0) is very sensitive to the choice of the 
coupling, especially for small renormalized strange quark masses, the mass functions at finite chemical potentials are 
not. As explained in the previous section, this is due to the medium modification of the coupling, which also leads 
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Figure 2: Mass functions at different values of three-momentum (see text) as function of the renormalized strange quark mass 
at a chemical potential of /i = 400 MeV for the coupling ai(k 2 ) (left) and ctu(k 2 ) (right). 



to significantly smaller mass values at the Fermi surface. Note furthermore that the values of the mass functions at 
the Fermi surface in the chirally broken phase and the CFL phase are very close to each other. This leads to the 
conclusion that the dynamics near the Fermi surface, where gapped and ungapped propagators strongly differ, are 
not directly relevant for the dynamical mass generation. As expected, the mass functions at the renormalization scale 
are comparable to the values in the chirally broken vacuum, which confirms that v is already sufficiently above the 
scale of dynamical mass generation and the Fermi energy. 



C. Dependence of the gap functions on the renormalized strange-quark current mass 

In this section we present results for the gap functions at the Fermi surface, i.e. at p4 = 0, and selected values of the 
three-momentum. The renormalization-point independent gap functions (|15f) for the triplet, doublet and anti-doublet 
channel (cf. Eq. 124(1 ) are evaluated at their corresponding Fermi momentum (see section lllBll . The functions 4>i^x0ud)i 
0i,l(«ds) and 0i,i(&s)> corresponding to M2, M4 and M3, given in the Appendix, are evaluated at (j>f,1i +Pf,l 2 )/2) 
(.Pf,1i + Pf,i 2 + 2p_F,i 3 )/4 and pf.i 3 , respectively. Due to this, the functions are considered at momenta corresponding 
to the pairing quasiparticles. This also allows to recover the results at the Fermi surface for the CFL phase in the 
chiral limit. The corresponding numerical results are shown in Fig. [21 for the coupling ai(k 2 ) and in Fig. 0] for the 
coupling au(k 2 ). Note again that due to the phase choice the gap functions 4>a, % and 4>d,i are real and 4>B,i and <pc,i 
imaginary at p^ = 0. 

All functions, apart from 4>i t i(bs), evolve towards the corresponding 2SC solution. However, 4>i.i^s) already shows 
that the transition must be first order. Furthermore, for the gap functions 4>i t i( u d) and ^,3, describing non-strange 
pairing, and <po vanish, as expected. Finally, we see for the other gap functions, relevant to strange pairing, that 
<pB is the most varying and 4>d is non- vanishing, which is in line with the interpretation of the gap functions given 
111 section UTO It again is obvious that the sensitivity of the gap functions on the used coupling is much weaker than 
those of the mass functions. 



D. Dependence of mass functions on the chemical potential 



We proceed by discussing the dependence on the chemical potential and treat the case of the mass functions first. 
Fig. [S] shows these mass functions in the doublet channel for two different renormalized strange-quark current masses, 
evaluated at the Fermi momentum and at the renormalization scale, respectively. 

For small enough chemical potentials, the 2SC phase is preferred and the doublet channel corresponds to decoupled 
strange quarks in the truncation used. Depending on the value of the renormalized strange quark mass, the strange 
quarks may not condense for small enough chemical potential, i.e. not develop a Fermi surface. This can be seen 
from Fig. [S] for a renormalized strange quark mass of m s = 200 MeV. In this case, we evaluate the mass function 
at vanishing momentum and the onset of strange quark condensation is reflected as a kink in the curves. With 
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Figure 3: Gap functions at selected values of the three-momentum for different channels (see text) as a function of the 
renormalized strange-quark mass at a chemical potential of fi = 400 MeV and for the coupling ai(k 2 ). 



rising chemical potential, the system undergoes a phase transition into the CFL phase. This effects the value of the 
mass function on the Fermi surface only slightly and is shown as a gap in the plots. For the coupling aji{k ) at 
m s = 100 MeV we find a direct transition of non-condensed strange quarks into the CFL phase. 



9 




"0 25 50 75 100 125 150 175 200 
m (v=2GeV) [MeV] 




u 25 50 75 100 125 150 175 200 
m (v=2GeV) [MeV] 




-100, 



25 50 75 100 125 150 175 200 
m (v=2GeV) [MeV] 



20 
15 

> 10 


-5 
-10, 





1 1 


1 1 1 1 1 




K 












— "i*B 














































1 1 1 1 1 1 1 



25 50 75 100 125 150 175 200 
m (v=2GeV) [MeV] 




"0 25 50 75 100 125 150 175 200 
m (v=2GeV) [MeV] 




"0 25 50 75 100 125 150 175 200 
m (v=2GeV) [MeV] 



Figure 4: Gap functions at selected values of the three-momentum for different channels (see text) as function of the renormalized 
strange-quark mass at a chemical potential of fi = 400 MeV and for the coupling au(k 2 ). 



The values of the mass functions at the renormalization scale only show a slight dependence on the chemical potential 
and are comparable to those in the chirally broken vacuum, which again reflects the fact that the renormalization 
scale is well above the dynamical chiral symmetry breaking scale and the chemical potential. On the other hand, 
the values of the mass functions at the Fermi surface already at a chemical potential of /i = 1 GeV are close to their 
values at the renormalization scale. As a result, dynamical chiral breaking is suppressed and the mass function is only 
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Figure 5: Mass functions on the Fermi surface and at a renormalization scale v = 2 GeV for fixed renormalized strange quark 
mass in the vacuum as function of the chemical potential for the coupling ai(k 2 ) (left) and au(k 2 ) (right). 



weakly dependent on the momentum below the renormalization scale. 



E. Dependence of gap functions on the chemical potential 



For completeness we also present results for the dependence of the gap functions on the chemical potential. In 
Fig. [S] we show the results for the gap functions in the CFL phase at a renormalized strange quark mass of m s (y — 
2 GeV) = 200 MeV and for the coupling aj(k 2 ). As described above, the 2SC phase is preferred for smaller chemical 
potentials and we find again a visible jump in the 4>i,i(ud) functions at a certain chemical potential. The gap functions 
in the CFL phase are remarkably insensitive in the chemical potential. Only the 4>b and cf>n functions evolve towards 
zero, which again reflects that the relevant values of the mass functions also become smaller. 



F. Pressure difference and critical strange-quark mass 



We now turn to the main result of this investigation: The determination of the critical value of the strange-quark 
current mass. Above this mass the Fermi surfaces are so far separated that pairing of up and down quarks with 
strange quarks is no longer energetically preferred. Based on the CJT-formalism we determine the pressure difference 
of CFL and 2SC using Eq. ©. The results as a function of the renormalized strange quark current mass for different 
chemical potentials and for the couplings employed is shown in Fig. [7\ One sees, as expected, that for small masses 
the CFL phase is preferred, and that there is a critical value of the strange-quark mass where the CFL phase becomes 
energetically disfavored. 

We like to emphasize that we are no longer able to find a solution for the CFL phase in case the 2SC phase becomes 
favored. This is considered as a consequence of the numerical method to solve the truncated DSE being the functional 
derivative of the truncated CJT action. It turns out that we always only find the global minimum of the CJT action 
as long as the local minimum is not protected by a higher symmetry. The latter is the case for the 2SC phase, if the 
CFL one is preferred. Nevertheless we can judge from the behavior of the gap functions that the transition is first 
order. 

In Fig. [S] the results for the critical value of the renormalized strange-quark mass as a function of the chemical 
potential are given and compared to the range of the physical strange-quark current mass as determined by the 
particle data group |l9|. (As stated above, the difference between this mass in the MOM and the MS regularization 
scheme is negligible compared to the experimental uncertainty.) As can be seen, the approach taken here predicts 
that the physical strange-quark current mass is very likely too small for allowing a 2SC phase at zero temperature for 
any chemical potential. This result is remarkably stable against the variation of the running coupling a s (k 2 ). 
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Figure 6: Gap functions on the Fermi surface for different channels (see text) at m a (y = 2GeV) = 200 MeV as a function of 
chemical potential for the coupling ai(k 2 ). 



IV. CONCLUSIONS AND OUTLOOK 



We have studied the quark propagator in the 2SC and CFL phase at zero temperature for different values of the 
strange-quark current mass in a fully self-consistent Dyson-Schwinger approach to color-superconductivity. Due to 
the medium modification of the interaction we find the 2SC phase to be disfavored at any physically relevant chemical 
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Figure 7: Pressure difference between 2SC and CFL phase as function of the renormalized strange quark mass at chemical 
potentials of n = 300 MeV, 400 MeV, 500 MeV for the coupling a/(fc 2 ) (left) and aii{k 2 ) (right). 
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Figure 8: Critical renormalized strange-quark mass as a function of the chemical potential for the couplings ai(k 2 ) and an (k 2 ) 
and the range stated by the particle data group (shaded band). 



potential for the physical value of the strange-quark current mass. This result is robust against variation of the 
running of the strong coupling in the infrared within the given uncertainties. Since the CFL phase is only mildly 
influenced by neutrality conditions we expect the result to hold even when neutrality is required. The critical value of 
the strange-quark current mass should then be similar or even slightly larger. This will be explored in further studies. 

Given the present results a further investigation of the CFL phase by the inclusion of the Meissner effect and 
Goldstone-boson contributions is of interest since both are expected to decrease the gap functions. Also the temper- 
ature dependence and therefore the structure of the QCD phase diagram at large densities is within the reach of the 
type of investigations presented here. 
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Appendix: THE COLOR-FLAVOR STRUCTURE OF THE CFL PHASE 



We generalize the ansatz for the CFL phase with two degenerate quarks as given in |4j by choosing the matrices 



Pi 



Mi 



( <5a + S i2 8 i2 5u 

Si5 5i5 $i3 



'hi 



Si7 



',8 



V 

/ Sn + 5 i2 S i2 S i4 , 

<h'5 &ih Sis 



Si7 



(In 



$i7 



\ 



Si8 
Si7 J 



which fulfill the requirements of Eq. 1221) . The basis is defined by 

{(r, u), (g, d), (b, s), (r, d), (g, u), (r, s), (6, u), (g, s), (b, d)}, 
with r, g, b denoting the color and u, d, s the flavor of the quarks. 
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